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ABSTRACT 

We present a stellar population synthesis study of a type I luminous infrared galaxy (LIRG): 
IRAS F13308+5946. It is a quasar with absolute magnitude M, = -22.56 and has a spectral feature of a Seyfert 
1.5 galaxy. Optical images show characteristics of later stages of a merger With the help of the stellar syn- 
thesis code STARLIGHT (Cid Fernandes et al. 2005) and both Calzetti et al. (2000) and Leitherer et al.'s (2002) 
extinction curves, we estimate the past infrared (IR) luminosities of the host galaxy and find it may have ex- 
perienced an ultraluminous infrared galaxy (ULIRG) phase for nearly 300 Myr, so this galaxy has probably 
experienced a type I ULIRG phase. Both nuclear starburst and active galactic nuclei (AGN) contribute to the 
present IR luminosity budget, and starburst contributes ~70%. The mass of supermassive black-hole (SMBH) 
is Mbh = 1.8 X 10^ Mq and the Eddington ratio Lboi/iEdd is 0.12, which both approximate to typical values of 
PG QSOs. These results indicate that IRAS Fl 3308+5946 is probably at the transitional phase from a type I 
ULIRG to a classical QSO. 
Subject headings: infrared: galaxies — galaxies: evolution — galaxies: starburst — 

galaxies: quasars: individual (IRAS Fl 3308+5946) — galaxies: stellar content 



1. INTRODUCTION 

Luminous infrared galaxies (LIRGs) and ultraluminous in- 
frared galaxies (ULIRGs) are galaxies with infrared (IR) lu- 
minosities Lir =L(8-1000^m) > 10" Lq andLiR > 10'^ L©, 
respectively. Their IR luminosities and space densities are 
both comparable to QSOs (Sanders & Mirabel 1996). These 
galaxies always show signs of tidal interaction and merger, 
and the interaction rate is observed to increase with IR lu- 
minosity (Zou et al. 1991; Clements et al. 1996; Murphy et 
al. 1996; Wu et al. 1998; Veilleux et al. 2002). Nearly all 
ULIRGs are strong interaction or merger systems (Kim et al. 
2002). 

It is widely accepted that galaxy interactions and mergers 
trigger extreme nuclear activity, as well as more widespread 
starburst (Toomre & Toomre 1972; Larson & Tinsley 1978). 
The molecular gas concentrations of ULIRGs are in their 
central kpc regions (Downes & Solomon 1998) and have 
the ability to form stars with densities comparable to those 
in elliptical galaxies (Tacconi et al. 2002). Kormendy & 
Sanders (1992) proposed that ULIRGs may evolve into el- 
liptical galaxies through merger induced, dissipative collapse. 
Structural, kinematic and photometric properties of ULIRGs 
show that they are originated from major mergers and they av- 
eragely fall on the fundamental plane of moderate-mass ellip- 
ticals (stellar mass ~ lO^'^-lO^M©), but are well offset from 
giant ellipticals, suggesting that ULIRG mergers are ellipti- 
cals in formation (Genzel et al. 2001; Tacconi et al. 2002; 
Dasyra et al. 2006a). 

Observations revealed that a large fraction of 
LIRGs/ULIRGs show spectral characteristics as Seyfert 
galaxies (Wu et al. 1998; Kewley et al. 2001). The fraction 
increases dramatically with IR luminosity, and among 
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ULIRGs, the fraction is about 50% (Kim et al. 1998a; 
Veilleux et al. 1999; Cao et al. 2006; Yuan et al. 2010; 
Nardini et al. 2010), or even higher (^70%; Nardini et al. 
2010). Though many of these sources are dominated by AGN 
in their bolometric luminosity (Boiler et al. 2002; Nandra & 
Iwasawa 2007), in most cases the predominance of starburst 
over AGN is proposed (Gu et al. 1997; Lutz et al. 1999; 
Franceschini et al. 2003), and Nardini et al. (2008) proposed 
a fraction of ^ 85%. 

From millimeter-wave CO observations and optical spectra 
analyse of 10 ULIRGs, Sanders et al. (1988) proposed that 
ULIRGs represent the dust-enshrouded stage of QSOs. They 
proposed a classical evolution scenario that two gas-rich spi- 
rals merge (or interact) first and the funneled gas toward the 
merger center triggers nuclear starburst before the ignition of 
a dust-enshrouded AGN. When dust has been consumed or 
swept away under the radiative pressure of AGN and super- 
novae, an optical quasar would appear Canalizo & Stockton 
(2000, 2001) selected a sample of 9 low-redshift QSOs which 
fall onto the intermediate position between the regions occu- 
pied by ULIRGs and QSOs in the far-infrared (FIR) color- 
color diagram. All these 9 transition QSOs are undergoing 
tidal interactions and 8 are major mergers. All of them also 
show strong recent star formation activity within 300 Myr 
They proposed a model involving a dust cocoon that initially 
surrounds the QSO nuclear regions, which can account for all 
the observed and derived properties in transition objects, and 
it is consistent with the idea of Sanders et al. (1988). They 
suggested that either at least some ULIRGs evolve to become 
classical QSOs, or some QSOs are born under the conditions 
as ULIRGs and their lifetimes as QSOs last < 300 Myr Most 
of these transition QSOs are also ULIRGs. As Kawakatu et al. 
(2006), we refer to QSOs and Seyfert 1 galaxies that selected 
from ULIRGs as type I ULIRGs (or IR QSOs in Zheng et al. 
2002). Kawakatu et al. (2006) and Hou et al. (2009) find that 
type I ULIRGs are the early phase of black-hole (BH) growth 
and they are QSOs in formation. 

The case may be as Colina et al. (2001) had proposed that 
high-luminosity QSOs would be the end point in the merging 
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process of massive (> L*) disk galaxies, while cool ULIRGs, 
which would be the end product in the merging of two or more 
low-mass (0.3L*-0.5L*) disk galaxies, would not evolve into 
QSOs. This scenario is supported by many studies (McLeod 
& Rieke 1994; Ambas et al. 2003; Dasyra et al. 2006a). 

In this paper we present the result of stellar population 
synthesis of IRAS Fl 3308+5946 and estimations of past IR 
luminosities. IRAS Fl 3308+5946 is a LIRG at present and 
shows a spectral characteristic of a Seyfert 1.5 galaxy, which 
together suggest a potential transition phase between type I 
ULIRGs and classical QSOs. It shows clear absorption lines 
and young stellar population features for spectral fitting. Ob- 
servations from IRAS, 2MASS, SDSS and FIRST are pre- 
sented in §2. We describe the process of extracting the star 
formation history (SFH), emission-line fitting, aperture cor- 
rection and the extinction curve in §3. In §4 we present some 
estimations like the current star formation rate (SFR), the BH 
mass, the Eddington ratio, etc., based on the emission-lines 
fitting. The past IR luminosities are estimated in §5. Conclu- 
sions and discussions are given in §6. We summarize in §7. 
Through out this paper we adopt a cosmology of //o=71 km 
S-' Mpc-i, nM=0.29 and r2A=0.71. 

2. OBSERVATIONS 

IRAS Fl 3308+5946 was observed by the Infrared As- 
tronomical Satellite {IRAS), Two Micron All Sky Survey 
(2MASS) All-Sky Extended Source Survey (recorded as 
2MASX J 13323783+5930538), Sloan Digital Sky Survey 
(SDSS, recorded as SDSS J133237.94+593053.7) and Very 
Large Array (VLA) Faint Images of the Radio Sky at Twenty- 
centimeters (FIRST) survey. It hasn't been explored by UV 
or X-ray observations by far 

2.1. Mid- and Far-Infrared 

The mid-to-far IR fluxes of IRAS Fl 3308+5946 are taken 
from the IRAS Faint Source Catalog (ESC) Version 2.0. Flux 
densities at three wavebands, 12, 25 and 100 /im, are upper 
limits. Only the one at 60 /im, which is 0.2615 Jy, has arecord 
with high quality (FQUAL=3). Therefore, the IR luminosity 
(LiR(8-1000/im)) is calculated by 



FiR(l-lOOOAim) = (1.75 ±0.25)FFiR(40-120Aim) (3) 



LiR(l-1000A*m) = 4TrDlFiR(l-l000^im)L(: 



(4) 



LiR(8-1000//m) « 2L, 



'60 /J. m ; 



(1) 



(Lawrence et al. 1989; Bushouse et al. 2002; Arribas et 
al. 2004; Wang et al. 2008), where Leo^m is the luminosity 
at 60 fim. The derived IR luminosity is LiR(8-1000/im) = 
lO'^'^^LoB For comparison, we also calculate IR luminosity 
using an IR spectral energy distribution (SED) model. Since 
IRAS Fl 3308+5946 is not only a starburst(SB) galaxy, it also 
has a type 1.5 AGN (see §3.4), we adopt the SED model es- 
tablished by Siebenmorgen et al. (2004) for galaxies with both 
AGN and starburst components. After the SED model was 
scaled to the observed 60/im flux, the integrated IR luminos- 
ity is calculated to be LiR(8-1000^m) = 10"^^ L©, exactly the 
same as the one above. 

We also calculate the IR luminosity in the range 1-1000 
fiva through 60 /im and 100 /im fluxes (Lonsdale-Persson & 
Helou 1987; Calzetti et al. 2000): 

FnR(40-120/im) = 1.26 x 10-'4{2.58/6o^n, + /iO()Mm}[W m'^] 

(2) 

^ The luminosity in tlie wavelengtli range 1-8 fim is estimated to be less 
than Ri5%-10% of the total and negligible (Calzetti et al. 2000). 



The derived IR luminosity LiR(l-1000/im) = lO'^^^o^^L© 
can be used as upper limit, and it is consistent with the result 
only adopting /eo/jm- Therefore, we adopt LiR(8-1000/im) = 
lO^'^^L© as the IR luminosity throughout this paper We 
don't do /T-corrections on IR fluxes, because the correction 
factor is small (~5% adopting SB SED model from Sieben- 
morgen et al. 2007) enough to neglect, compared with the un- 
certainty of IR luminosity calculation and the scatter among 
different /^-correction methods. 

2.2. Optical and Near-infrared 

IRAS Fl 3308+5946 was also observed by SDSS imaging 
camera and spectrograph. The u-, g-, r-, i-, z-band images and 
the spectrum are from the SDSS Data Release 5 (DR5). Each 
band was observed with exposure time of 53.9 s, and the pixel 
size is 0.396" on the sky. The spectrum has been corrected 
for Galactic reddening of ^(B-y) = 0.015 and redshift of z = 
0.171, so the rest-frame wavelength coverage is 3258-7852A. 

The SDSS spectrum shows a typical feature of a type I 
AGN. The /-band absolute magnitude of M,- = -22.564 sat- 
isfies the quasar criterion of M, < -22.0 defined by Schnei- 
der et al. (2007). However, the B-band absolute magnitude 
is Mb = -21.1 (transferred from Mb = -21.8 which adopts 
Hq = 50 km s"' Mpc"', Veron-Cetty & Veron 2006), and 
it doesn't satisfy the criterion for quasars of Mb < -22.2 
(Ho = 71 km s-i Mpc-i) defined by Schmidt & Green (1983). 
The larger extinction at B-band could explain the inconsis- 
tency between the two band (see also §3.6 ). 

There is obvious deviation between SDSS and IRAS posi- 
tions. The IRAS 2.45-cr uncertainty ellipse (corresponding to 
95% probability enclosure) with semi-major and semi-minor 
axises of 21" and 7" does not embrace the optical (SDSS) po- 
sition (see Fig.[T]i. The optical source only locates within the 
5-0- uncertainty ellipse of IRAS position. Since there is not 
any other bright sources within IRAS uncertainty ellipse in 
either optical and near-infrared (2MASS) bands, and IR phe- 
nomenons are most related with galactic activities, we identify 
the object detected by IRAS, 2MASS and SDSS as the same 
one. In §2.3, the well-known radio-FIR correlation will con- 
firm such identification. Hereafter, we refer to the object as 
IRAS Fl 3308+5946. 

Fig. |2] (left) gives the M-band image and the 3" diameter 
fiber of the spectrograph, which covers ^10 kpc of the galaxy. 
Two color-map images of u-g and g-r are also presented. 
The field of view (FOV) of these images are all 25"x25". 
These images are made after world coordinate system match- 
ing. Galactic reddening correction and /T-correction (IDL rou- 
tine from Blanton et al. 2003, v4_l_4). Because the seeings 
at U-, g-, r-bands are 1.3", 1.3" and 1.1" respectively, we have 
to match the seeings by convolving the r-band image with a 2- 
D Gaussian kernel with a = 0.29" before making g-r image. 
The g-r image shows a color gradient that the central region 
of galaxy is quite blue. It could be explained by circumnu- 
clear starburst and/or AGN activity. After converting SDSS 
g-r color to B-V color (Jester et al. 2005), we obtained the 
B-V color varies from ^0.25 at the center, to ^0.5 at the edge 
of the fiber aperture and to > 1 at the outskirt of the galaxy. If 
these variations are due to different stellar populations, the 
ages of the stellar populations may range from ^300 Myr to 
~800 (assuming solar metallicity) inside the aperture, and to 
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wavelength. It models M\ by a combination that 
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Fig. 1 . — Source identification from IRAS and SDSS. Tlie red ellipse 
is IRAS uncertainty ellipse with 2.45-cr confidence for each axis, which 
is plotted on SDSS r-band image. Blue circles represent SDSS galaxies 
and magenta ones SDSS stars. The cyan circle marks the FIRST source 
in the field of view, and coincides with the galaxy IRAS F13308+5946. 



>IOGyr at the galaxy outskirt, adopting the stellar popula- 
tions evolution model of Bruzual & Chariot (2003; hereafter 
BC03). Considering local dust extinction, the g-r color may 
be reddened, so the ages obtained above should be somewhat 
younger Generally, the u-g color distribution is consistent 
with g-r color according to BC03. However, unlike g-r color, 
the bluest region on the u-g image deviates from the galaxy 
center. It could be attributed to the merger process. 

From Fig. [T] and Fig. |2] we can find characteristics of later 
stages of a merger. Adopting 2MASS standard isophotal pho- 
tometry, the //-band absolute magnitude is Mh =-25.7, which 
is 4 times as luminous as L* galaxies (Mh = -24.2; Colina et 
al. 2001). Therefore, the galaxy could be a merger system of 
massive (> L*) galaxies. 

2.3. Radio 

The only one FIRST source in the FOV of Fig.[T]is marked 
with a cyan circle. The radio source coincides with the opti- 
cal quasar with a spacial separation of 0.63". The 1.4 GHz 
continuum integrated flux density detected by FIRST is 0.93 
mJy. It follows the well-known radio-FIR correlation which is 
confirmed by Yun et al. (2001) by the overwhelming majority 
(>98%) of 1809 IR-selected galaxies (^fio^m > 2Jy). More- 
over, the derived parameter log/?* = log ■^'^°°"'!f ""' = 0.423 

(Sramek & Weedman 1980) is much lower than 1, indicating 
that it belongs to radio-quiet populations (Stocke et al. 1992). 
The coincidence between optical and radio position and its 
radio-FIR relation support that the optical and the IR source 
is the same one. 



3. SPECTRAL ANALYSIS 

3.1. Spectral Synthesis with STARLIGHT 

We use the spectral synthesis code STARLIGHT (Cid Fer- 
nandes et al. 2005) to derive stellar populations of the host 
galaxy. STARLIGHT fits the observed spectrum 0\ with a 
model spectrum M\ which is made up of a pre-defined set 
of base spectra. There is no need to give the value of pa- 
rameters initial guess. It carries out the fitting with a simu- 
lated annealing plus Metropolis scheme to yield the minimum 



X 
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= X)a[('^^~-'^^)^^]^' where w^ is the error in 0\ at each 



where 7;,a = b\j ® G(v*, ai,), b\j = (j^^) is the normahzed 

flux of the /* spectrum, B\j is the /'' component of base 
spectrum, Bx^j is the value of the /'' base spectrum at the 
normalization wavelength Ao, G(vv, cr*) is the Gaussian distri- 
bution centered at velocity v^, and with dispersion a^, of the 
line-of-sight stellar velocity, Xj is the fraction of light due to 
component^ at Ao, r\ = 10""'*^'*^"'*''^ is the global extinction 
term and represented by Ay- The nebular and AGN emission 
lines can be obtained by subtracting the model spectrum from 
the observed one as E\ = 0\ -M\. 

STARLIGHT uses random Markov Chains which needs to be 
appointed an integer seed to generate random numbers during 
the fitting, and if it uses different seeds, it will give different 
results even if the parameters are configured identically. Such 
difference will not change the overall populations distribution 
significantly, but individual component .jcj may vary ^10 per- 
centage in some cases. To obtain a statistically reliable result, 
we carried out fittings with a set of seeds generated by Monte 
Carlo sampling and adopt the mean value over all seeds for 
every parameter 

3.2. Parameters Determination 

We take simple stellar populations (SSPs) from BC03 and 
a power-law spectrum of AGN (if needed) /a oc A"^ as our 
base spectra. To minimize the parameter space and to make 
the result more reliable, we adopt spectral templates with 10 
ages (0.005, 0.025, 0.1, 0.29, 0.5, 0.9, 1.4, 2.5, 4 and lOGyr) 
computed with "Padova 1994" evolutionary tracks (Alongi et 
al. 1993; Bressan et al. 1993; Fagotto et al. 1994a,b; Girardi 
et al. 1996) and Chabrier (2003) initial mass function (IMF). 

As a type I galaxy, the spectrum of IRAS Fl 3308+5946 is 
possible to contain a power-law component, though it is not 
obvious in the observed spectrum. Since STARLIGHT can in- 
clude a power-law spectrum in the template base, we can carry 
out test fittings to confirm whether it exists and find out the 
value of spectral index a a ■ A widely adopted AGN power-law 
slope (Fl, ex J/"") over the optical and UV region is a^, = -0.5 
(Richstone & Schmidt 1980), while other studies had given 
the values from -1 to (Natali et al. 1998, and references 
therein), corresponding to -2 < aA < -1- We carry out test 
fittings with aA varies from -3.0 to -0.1 at intervals of 0.1 to 
search for the power-law index, and additional non-power-law 
fittings with a\ = 0. 

We optimize the fitting with single metallicity SSPs. To 
search for the best-fit metallicity, we carry out metallicity test 
with the six metalHcities (Z = 0.0001, 0.0004, 0.004, 0.008, 
0.02, 0.05) of BC03 model for every power-law index. After 
that, we adopt the best-fit metallicity and power-law index for 
the formal fitting. 

Fig.[3]shows the test fitting results with different power-law 
indices a\ and the six metalHcities. We evaluate the fitting 
quality by the minimum x^, which is averaged over 25 fittings 
with different seeds (randomly and uniformly distributed be- 
tween -100000 and 100000). In Fig. [2 the minimum x^ set 
reaches its bottom at Z=0.008 and a.\ = -2.0. This is con- 
sistent with other studies that metalHcities near solar value 
prevail in starburst galaxies (Tadhunter et al., 2005; Pellerin 
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Fig. 2. — SDSS K-band image and u-g, g—r color maps of IRAS Fl 3308+5946. The red circle represents coverage of the 3" diameter fiber aperture on the host 
galaxy. The FOV of the three images are all 25" X 25" . 
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Fig. 3. — Power-law index and metallicity test. Every test fitting adopts 
one of the six metallicities, a certain power-law index and 10 SSPs, and is 
evaluated by averaged x's over 25 individual fittings with different seeds. 
The left panel shows the averaged x's versus metallicities (Z) and power-law 
indices (a). The right panel gives the color code of x^- 

& Robert 2007), and also that the power-law continuum slope 
falls around a\ = -2 (Francis 1996; Vanden Berk et al. 2001; 
Letawe et al. 2007). Moreover, the power-law index doesn't 
affect the synthesis very much. We learn from Fig. |3] that 
if we fix the metallicity, the variation of x^ among different 
power-law indices is not significant. Its influence on power- 
law component fraction is showed in Fig.|4] These data points 
are all from the test fits above, but only with subsolar metal- 
licity (Z = QAZq). The red unfilled diamonds with a fixed 
power-law index represent the AGN component fraction at 
4020A of 25 different integer seeds. We can see from the 
figure that in the normal range of power-law index in liter- 
ature, say, -2 < a>, < -1, the average AGN fraction keeps 
almost constant, and the dispersion of the average AGN frac- 
tion is far less than the one from different seeds. This result 
indicates that the form of AGN component doesn't affect the 
synthesis very much, so neither it affect the stellar component 
significantly. 

Therefore, we take the parameters ax = -2.0, Z = 0.008, 
and 10 SSP ages for formal fitting. 

3.3. Formal Fitting and Stellar Populations 

The formal fitting involves 100 independent fittings with 
different seeds. Stellar absorption lines are given three times 
the weight larger than the continuum to emphasize the de- 
tail of stellar component. The recognized emission lines like 
[O ii] AA3726H-3729 ([O u] A3727), Ne m A3869, H/3h-[0 m] 




Fig. 4. — AGN fractions in test fittings among different power-law indices. 
The red unfilled diamonds with a fixed power-law index represent data points 
of 25 different seeds. The blue filled circles represent the average AGN frac- 
tion of the 25 points. 

AA4959H-5007, HaH-[N n] AA6548H-6583, etc., together with 
Nal D5890 ISM absorption line are all masked before fitting. 
The blended Fe u emission-line regions are also excluded, 
which are assigned to be 4270^700 A and 5100-5600 A. 
These masked regions are illustrated in Fig.|5]represented by 
the gaps of the blue curve on the left-top panel. The blue curve 
is the error of the observed spectrum, and the inverse of error 
is the weight at each wavelength in spectral fitting. The ob- 
served spectrum Ox and all SSP templates are normalized at 
wavelength Aq = 4020 A before fitting. 

The outcome fraction for each stellar component varies 
within ~10 percentage point. Such variation is insignificant 
for main populations. We take the mean value for each pa- 
rameter as the formal fitting result, such as Xj, v*, (Ji,, y^, etc. 
These parameters are all summarized in Fig. |5] The left-top 
panel shows the observed spectrum Ox {green), the model Ma 
(red) and the error {blue). The left-bottom panel gives the 
residual spectrum Ex = Ox-Mx- Light- weighted stellar pop- 
ulation fractions Xj are shown on the right-top panel, on which 
the left black bar represents the power-law fraction without 
a concrete age. Mass-weighted population fractions /ii are 
shown on the right-bottom panel. We learn from Fig. |5] that 
the stellar populations of the galaxy consist of two main pop- 
ulations: a younger one <500Myr (whose 500 Myr popula- 
tion is insignificant) and a mediate-to-old > 1 Gyr (we define 
young stellar populations as those with ages < 300 Myr, old 
populations with ages > 1 Gyr, and intermediate-age popula- 
tions between them). Follow Cid Fernandes et al. (2005), we 
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Fig. 5. — Spectral synthesis of IRAS F13308+5946. Left-top: The observed spectrum 0\(green), the model spectrum Mx (red) and the error 
spectrum (blue) with the gaps meaning the masked regime and the thrice weighted absorption lines, left-bottom: The residual spectrum Ex 
(purple). Right: Light (top) and mass (bottom) weighted stellar population fractions Xj and ^j, respectively. The inserted panel on the right 

marks the ages of the stellar population templates. The flux intensities of the left two panels are both normalized at 4020A by 4.549872 x 

10"'*(ergs s"' cm"-). 



calculated the mean ages weighted by light, 



(logf^)L = "^Xjlogtj 



and by stellar mass 



M. 



(log4)M = ^Ai>logf/ 



(6) 



(7) 



The derived mean ages are (logf*)^ = 8.2 and (logf^)^ = 9.7. 
These figures demonstrate again that the galaxy mass is dom- 
inated by old populations with ages of ^ 10'*^ years, and a 
latest starburst took place about 10** years ago whose massive 
stars contribute a large fraction of luminosity. The mean ve- 
locity dispersion is a^, = 223 km s~', which is typical of giant 
ellipticals. 

3.4. Emission-Line Fitting 

The residual spectrum Ex (the purple spectrum on the left- 
bottom panel) in Fig. |5] can be used to measure emission 
lines, since the stellar component and the AGN continuum 
have been removed. Fe n emission regions are not included 
in emission-line fitting. Emission lines are modeled by the 
SPECFIT task in the IRAF-STSDAS package. 13 compo- 
nents are fitted simultaneously, including the single emission 



line [O ii], [N n] and [S n] doublets, narrow and broad compo- 
nents of H/3, Ha, [O III] A4959 and [O m] A5007 (represented 
by H/3n, H/3b, and so on). Each line as well as the narrow and 
broad components are fitted by a single Gaussian profile. The 
flux ratios of [O \i\\n, [O \i\\b and [N n] doublets are fixed at 
their theoretical values. The relative positions of [O hi], [N 
ii] and [S n] doublets are constrained by their laboratory val- 
ues. FWHMs are constrained to be same for each doublets 
like [O ui]a', [O ni]^, [N u], [S n] and also H/3, Ha narrow 
and broad components. The emission-line properties without 
aperture or extinction correction are displayed in Table 1 . The 
uncertainties in columns (2) and (3) are given by the SPECFIT 
task. Fitting results are shown in two wavelength ranges in 
Fig. |6] The narrow components of Ha and H/3 have FWHMs 
<1000km s"\ and their broad components are >7000km s~'. 
It can be classified as a quasar with spectral type of Seyfert 
1 .5 galaxy, for both narrow and broad components of the two 
Balmer lines are significant. 

Greene & Ho (2005) found a tight correlation between 
Balmer emission-line luminosities and AGN optical contin- 
uum luminosity (little host galaxy contamination) at 5 100 A 
(L5100 = ALa at a = 5 100 A). We use this correlation to ex- 
amine whether our spectral decomposition is successful. We 
measure Lgioo from the power-law spectrum which is obtained 
from spectral synthesis. The model AGN spectrum gives 
^5100 = 10'*'""^ L0. The observed Hae and H/3b luminosi- 
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ties are Lhqb = 10 L© and Lh,3b = 10 L©. Using Greene 
& Ho's correlations, L5100 is calculated to be lO^'^^L© and 
lO'^^L©, respectively. They are both consistent with the 
5 100 A luminosity we modeled, so the spectral decomposi- 
tion of the AGN and host galaxy is reasonable. 



TABLE 1 
The emission-line properties of IRAS F 1 3308+5946. Fluxes are 

MEASURED WITHIN THE FIBER APERTURE, WITHOUT EXTINCTION 
CORRECTION, AND IN UNIT OF 4.549872 X 10""'(ergs S"' cm"^). 



Line Identification 


Flux 


FWHM(kms-') 


[U„JA3727 


3.96±0.23 


725.85±4L33 


H/3n 


L96±0.23 


587.42±35.85 


H/3b 


12.17±0.53 


7403.59±97.26 


Hon 


12.01±0.66 


587.42±35.85 


Has 


55.57±L06 


7403.59±97.26 


[Oi„]A4959n 


4.05±0.16 


779.80±30.76 


[0,„]A4959b 


2.83±0.17 


2879.82±184.07 


[Oiii]A5007n 


12.26±0.49 


779.80±30.76 


[Oiii]A5007b 


8.57±0.51 


2879.82±184.07 


[Nn]A6548 


4.29±0.21 


543.60±27.18 


[Nii]A6583 


12.69±0.62 


543.60±27.18 


[Sii]A6716 


2.41±0.23 


516.52±49.54 


[Sii]A6731 


L91±0.20 


516.52±49.54 



3.5. Aperture Correction 

We have showed in Fig. |2] that the observed spectrum is 
the flux from the 3" diameter fiber of SDSS spectrograph. 
Precise calculations must include the flux come from the 
whole galaxy. Table 2 lists the Galactic reddening corrected 
fiber magnitude (fiberMag, measured within the aperture of a 
fiber), Petrosian magnitude {petroMag, measured by the mod- 
ified form of the Petrosian system) and A'-correction factor for 
each band. The aperture effect is significant, for the difference 
between fiberMag and petroMag can be larger than 1 magni- 
tude at some bands. Since there's no lager aperture spectrum 
of the galaxy available by now, aperture correction is neces- 
sary in this work. A rough estimation of the aperture effect 
can be derived via 



A - Lpetro/i'fiber - 10 



-0.4(HiPM„,-mr,bcr) 



(8) 



(Hopkins et al. 2003). The correction factors for «-, g-, r-, /-, 
z-bands are 1.68, 2.2, 2.58, 2.60 and 2.57, respectively. At 
M-band, the aperture effect is the least, but it is still a factor 
of A = 1.68. From the spectral synthesis in §3.3, it is known 
that the luminosity of the galaxy is dominated by young stellar 
populations and the AGN, which both can be better traced by 
M-band rather than other four bands, so we carry out aperture 
correction using w-band's factor. 



TABLE 2 
Magnitudes from SDSS fiber and Petrosian photometry. 



mag 



u 



1 



fiberMag 
petroMag 
k-correction 



18.46 
17.90 
0.028 



17.86 
17.00 
0.166 



17.37 
16.34 
0.109 



17.06 16.84 
16.02 15.82 
0.0213 0.086 



Equation ([8]l is a simple aperture correction, because it 
doesn't distinguish different origins of the observed light. The 
3" fiber is of two orders larger than typical AGN scale (for 
narrow-line region, ^100 pc), so the 3" fiber and the Petrosian 
radius may both embrace the whole AGN region. However, 



the stellar content varies in the two apertures (we see Petrosian 
radius a larger aperture). If we want to derive the flux of a star- 
burst powered emission-line coming from the whole galaxy 
and adopt equation (O, we will amplify the AGN component 
and therefore underestimate the contribution from stars. 

The spectral synthesis method provides us an opportunity 
to improve the aperture correction method above. Having de- 
composed the observed spectrum into AGN power-law and 
stellar component, we can convert their fluxes at a given band 
to apparent magnitudes. In combination with njfiber and mpen-o, 
we give a new method for purely stellar component aperture 
correction. The derived formula is: 



A = 



-'*fiber 



JQ-0.4(mpa„-mAGN) _ J 

JQ-0.4(mr,bcr-«lAGN) — 1 



(9) 



where L*petro and L^gbei are fluxes of stars at a given band 
within Petrosian and fiber aperture, respectively; wagn is the 
apparent magnitude of AGN at a given band converted from 
the flux density of the power-law spectrum. 

3.6. Extinction Curve 

Large extinction happens in ULIRGs/LIRGs and re-emits 
UV-to-optical emission to FIR. The proper dust extinction law 
we use here is the one given by Calzetti et al. (1994; 1997; 
2000) and Leitherer et al. (2002). Calzetti curve is only ap- 
plicable to wavelengths longer than 1200A. Leitherer et al. 
(2002) extended the curve to 970A. At even shorter wave- 
lengths, extinction curve has larger systematic uncertainties, 
but we extrapolate it a little to 912A to cover the Lyman series 
limit. Their curves are identical at 1500A, and the differences 
at shortward wavelengths are minor. We adopt Leitherer et 
al.'s curve for 912A< A < 1800A, and adopt Calzetti curve 
for 1800A< A < 9000A. 

The intrinsic emission /^(A) can be recovered through the 
extinction curve k (A) via 



Fi(A) = Fo(A)10' 



f.)A[E,(B-V)k (\)-E,{BA/)k^] 



(10) 



where F„{\) is the observed spectrum, Es{B-V) is the 
color excess for the stellar and AGN continuum spectrum, 
lQQA[EAB^)k {\)-EAB^)ky] coiTesponds to the global extinction 
term r\ in equation (|5]l. EsiB-V) is calculated from Ay 
obtained from the spectral synthesis (Ay = 0.62) and it is 
Es{B-V) = 0.15. The color excess for nebular gas emission- 
lines is denoted by Eg(B-V) and is directly estimated through 
Balmer decrement, the line ratio Hq;n/H/3n = 6.13, whose in- 
trinsic flux ratio is (HaN/H/3N)o = 2.87 assuming temperature 
T = 10000 K and case B recombination (Osterbrock 1989). 
Eg{B-V) is derived to be Eg{B-V) = 0.71. 

When applying the Calzetti curve to B-band, the extinc- 
tion correction factor is ^^4, corresponding to 1.5 magnitude 
brighter, such the corrected absolute magnitude is Mb = -22.6, 
which satisfies the quasar criterion Mb < -22.2. 

4. SFR AND BLACK-HOLE MASS 

We can compare the observed IR luminosity with the cur- 
rent SFR. Because the Ha line contains emission from the 
AGN, we use [O u] A3727 to estimated the SFR. We take 
color excess Eg{B-V) = 0.71 for gas reddening and Calzetti 
curve for spectral line extinction correction. The current SFR 
is calculated by 

SFR([0 n])(M0 yr"') = (6.58 ± 1 .65) 
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Fig. 6. — Fits to the emission lines. Flux intensities are inlierited from tlie residual spectrum in Fig. [5] Tlie observed line profile (green-solid), the model profile 
{red-solid), each line and different components (black-dashed) are plotted in each panel, left: H/3+[0 hi] region; right: Ha+[N n]+[S u] region. 



X 10-42 L[o„](ergss-i) (11) 

(Kewley et al. 2004), where L[o „] is the extinction corrected 
(a factor of 45.6) luminosity of [O n] A3727. The derived SFR 
is 43. 8± 1 1 .OMq yr"' . We don't take aperture correction here, 
because we have learnt that the starburst takes place mainly 
at the center region of the galaxy (see §2.2), where the fiber 
aperture covers. Considering the star- forming region can be 
also found outside the aperture, we take such SFR as a lower- 
limit and the M-band aperture corrected SFR as an upper-limit, 
which is 73.5±18.4Mq yr'^ For starburst galaxies, SFR can 
be also estimated through IR luminosity (Kennicutt 1998, and 
references therein). As for IRAS Fl 3308+5946, since it con- 
tains both starburst and AGN components, the total IR lumi- 
nosity 10" '^^L© only suggests an upper limit of SFR, which 
is 62.9 Mq yx~^ . The spectral synthesis has given the starburst 
fraction of the galaxy, so the SFR can be estimated only by 
the IR luminosity from the starburst. Such calculation will be 
given in §5. 
The BH virial mass is calculated by means of Hq;b via 



Mbh = (2.0![1;^)x10''' 



-^Hae 



/FWHMea 
V lO^km s-i 



\ lO^^^ergs s ' 

2.06±0.06 

Mq 



0.55±0.02 



(12) 



(Greene & Ho 2005). Where Lhqb i^ the observed Has lu- 
minosity. We don't carry out aperture correction on Lhqb' 
because the 3" fiber aperture covers ^^10 kpc, including 
the broad-line region definitely. The derived BH mass is 
Mbh = 1 .8 X 10** M0. The Eddington ratio Lboi/i'Edd is calcu- 
lated adopting the bolometric luminosity Lboi ~ 9AL;j(5 lOOA) 
(Kaspi et al. 2000), where La(5 lOOA) is the AGN optical con- 
tinuum luminosity at 5100 A. The inferred Eddington ratio is 
Lboi/i'Edd = 0.12, which is typical of PG QSOs (median value 
0.24; Hao et al. 2005, hereafter Hao05) and far less than the 
one of IR QSOs (1.73) and NLSls (1.27). These results in- 
dicate the SMBH is at the late stage of growth. If the SMBH 
grew from a single seed, the duration of the growth can be 



estimated by the e-folding timescale, 

Mbh 



'acc — 



M 



:4X 10 



'Ory 



yr 



(13) 



(Haiman & Loeb 2001), where e = Loa\/Mc~ is the mass-to- 
radiation conversion efficiency and 77 = Lboi/^Edd is the Ed- 
dington ratio. For efficient energy production, e '^ 0.1 and 
the e-folding timescale is 4cc = 3.5 x 10^ yr. It is close to the 
duration of the starburst phase. This means that before the 
starburst activity took place, the BH mass was less by a factor 
of ^2.7, which was still massive enough to support a broad- 
line region. Therefore, the galaxy was probable to harbor a 
type I AGN during the starburst phase. 

5. FIR LUMINOSITY IN THE R\ST 

We have little information to determine the origin of the 
BH, that is, whether it grew from a single seed or from BHs 
merger when galaxies merged. Therefore, the AGN emission 
in the past is hard to estimate. Comparing with the uncertainty 
of AGN, we have better defined stellar evolution models (Ku- 
rucz 1992; Lejeune et al. 1997, 1998; BC03). Combining 
with the star formation history, we are able to trace back to 
the past and calculate the luminosity of the stellar component 
at a given epoch. Then, we could estimate the IR luminosity 
by calculating the UV-to-optical flux absorbed by dust. Be- 
cause there is no UV observation data available for this galaxy 
and the shortward end of the rest-frame spectrum is only lim- 
ited to 3258 A, we use model spectrum of UV-band for such 
calculation. 

To test the feasibility of our method, we firstly calculate the 
present FIR luminosity by reconstructing the UV-to-optical 
spectrum. Because BC03 SSP templates are in units of so- 
lar luminosity per angstrom per solar mass, we can accurately 
convert the present stellar mass Mcor_tot to a spectrum by 
means of the mass weighted fraction /ij. We reconstruct the 
model spectrum Fo{Xfy) covering the whole UV-to-optical 
wavelengths in our calculation (from 912A to 9000 A) through 

N 

Fo(A,fo)= [Mcor_tot^AijBA,j,to+ivXA)]10-''-4<^^-^''> (14) 
j=i 
where /q is the present time, Mcor_tot is the present stellar 
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mass obtained from the spectral synthesis and is derived to be 
5.5 X IO'^Mq after aperture correction through equation (|9]l 
(doesn't need extinction correction), Bajio are BC03 SSP 
templates without normalization, Fp{X) is a double power-law 
spectrum of AGN. The spectral indices of power-law we use 
here are given by Hatziminaoglou et al. (2008) that a = -l for 
A < 1250A and a = -2 for A > 1250A. Fig. 13 shows the re- 
constructed model spectrum {red-solid) superimposed by the 
observed one (green-solid). This model spectrum is used as 
F„(X) in equation ( fTOl ). 




\ (A) 

Fig. 7. — The reconstructed model spectram (red-solid) of the present state. 
The observed spectrum (green-solid) is superimposed on it. The blue-solid 
hne represents the host galaxy starlight, and the black-dashed line represents 
the double power-law spectrum. 



The light in UV-to-optical (912-9000A) absorbed by dust 
and re-emitted to IR is approximately calculated by 



.^IR,/ ■ 



9000 



912 



[Fi(X,t)-F,(X,t)]dX 



(15) 



where Fi(X,t) and Fo{X,t) have similar meaning as Fi(X) and 
Fo(X) in equation ( fTOl i. but they represent the spectra at look- 
back time t. After M-band aperture correction, the derived 
IR luminosity is Lirq = lO'^^^L©. It approximates to the 
observed Lir = lO'^^^L©. It demonstrates that our method 
is feasible in recovering the present FIR luminosity, so we 
can apply the method to the past. Meanwhile, it is conve- 
nient to calculate the contribution fractions from starburst and 
AGN. We only need to calculate IR luminosities from the two 
components separately. The starburst and AGN contribute 
istar.iR = 10"'^^ L0 andLAGN.iR = IO'^'^'^Lq, respectively. The 
starburst contributes ^ 70% of the FIR luminosity, nearly two 
thirds of the total FIR budget. Anyway, two components are 
both important. 

The current SFR can be estimated through the relation be- 
tween LFiR(8-1000^m) and SFR (Kennicutt 1998). Since the 
host galaxy contributes ~ 70% of the FIR luminosity, we ob- 
tain SFR = 43.3M0 yr"' . It is consistent with the one from [O 
ii] luminosity inside the aperture. Therefore, the derived stel- 
lar populations and AGN spectra, the dust extinction and the 
observed emission line fluxes are concordant with each other. 

The estimation of the past IR luminosity is carried out in a 
similar way. We just need to supplement the stellar mass loss 
of each population at a given look-back time t. Because the 
AGN emission in the past is unknown, we just estimate the 



stellar component contribution. The stellar spectrum (i.e. the 
host galaxy) in the past without extinction is produced through 

N 

F,(X,t) = Mcor_tot^ -p-U,uBx,u (16) 

where /*]_,„ is the present (fo) fraction of the remaining stel- 
lar mass to the initial mass of population j, /^j , is such 
fraction at a given time t, f-^to and /+.j., are in the range 
< /*.j.'oi./*j.' ^ 1- From Fig. |5](the /7offom-r/g/if panel) Since 
the recent starburst happened in the past 500 Myr, we estimate 
the IR luminosity of the past 25, 100, 290 and 500 Myr sepa- 
rately. 

Take 25 Myr ago as an example. When the 25 old popula- 
tion was just newly born, it was assigned an age of 1 Myr and 
the parameter /*.j., is equal to 1. The 100 Myr old population 
was 75 old then with /* j , = 0.7226. Other populations are also 
configured in the same way. Here, the reconstructed spectrum 
through equation ( fTSI l is used as the intrinsic spectrum Fi(X) 
in equation ( fTOl ). The extinction term A a is assigned to be as 
the present. The attenuated spectrum through Fi(A,f)10~"'*'^^ 
is used as /\,(A) in equation (fTOl ). Because dust extinction in- 
creases from optical to UV, we adopt the aperture correction at 
M-band derived from equation (|9]l. Calculated through equa- 
tions ([Tol l and ( fTSl l. the IR luminosity transferred from 912- 
9000A to IR at 25 Myr ago was Lir.25m = lO'^'^L©, which is 
ULIRG luminosity. If we reckon in the AGN's contribution, 
the IR luminosity should be higher. When applying to 100, 
290 and 500 Myr ago, the estimated IR luminosities from the 
host galaxy were Lir.ioom = lO'^^'^Lo, Lir,290m = lO'^^'L© 
and L1R.500M = 10"'*^Lo. These IR luminosities are shown as 
a history in Fig. [8] As has been demonstrated in section §4, 
the galaxy was probable to harbor a type I AGN during the 
starburst. Therefore, a type I ULIRG might appeared during 
the starburst phase and lasted for '--^SOO Myr. 




t (Myr) 

Fig. 8. — The probable IR luminosities during the stai'burst phase. The 
ULIRG phase may lasted for ~300Myr. 



6. DISCUSSION 
6.1. Fe II pseudo-continuum 

Many IR QSOs are extremely strong Fe n emitters (Zheng 
et al. 2002; Lipari et al. 2003). In IRAS F13308H-5946 
spectrum, we also find Fe n pseudo-continuum emission 
in wavelength region 4270^700A and 5100-5600A. We 
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have masked these two wavelength regions before spec- 
tral synthesis. Actually, we have attempted to subtract Fe 
n lines from 0\ before fitting. However, as the case of 
IRAS Zl 1598-01 12, IRAS F02065+4705, etc. (Zheng et al. 
2002), Fe u multiplets 37, 38 (4500-4680A) are relatively 
stronger than the Fe u multiplets 48, 49 (5 100-5400 A) com- 
pared with Boroson & Green (1992, hereafter BG92) Fe u 
template, and also the Veron-Cetty et al. (2004) template. 
Thus, the Fe u multiplets 37, 38 are left when the multiplets 
48, 49 are removed. We have also attempted to remove Fe 
n lines from the residual spectrum, but the result was not im- 
proved. Furthermore, including Fe n template as a component 
to fit 0\ does not bring any improvement at all. Therefore, 
we masked the Fe u pseudo-continuum regions before spec- 
tral synthesis, as have mentioned in §3.4. 

The Fe u A4570 (4434-4684A) flux is measured during the 
procedures above and the line ratio Fe n A4570/H/3 is 0.81. 
This ratio is typical of optical selected QSOs (BG92). We 
find that the Fe n continuum extension to [O n] A3727, H/3, 
etc. lines is hard to detect, so we measure these lines fluxes 
directly. 

6.2. Extinction Curve 

In this paper we don't adopt Cardelli et al.'s (1989) extinc- 
tion curve, because it is for cases which are similar to Milky 
Way and not applicable to starburst galaxies. We use the 
Calzetti curve and Leitherer et al.'s (2002) extension form, be- 
cause they are derived from star-forming regions and starburst 
galaxies. 

Dust in different components and stellar populations vary 
in temperature, ingredient, amount, geometry, etc., so differ- 
ent attenuation for each component should be adopted if we 
could distinguish them. However, as the spectra from the host 
galaxy and central AGN are coupled, and also there is only 
one extinction term in STARLIGHT, the extinction obtained 
from spectral synthesis is a combination value from both re- 
gions. Also, there may exist some optically thick star forma- 
tion regions, which may not be detected even by optical or 
near-infrared, but they still contribute to IR flux. In this case, 
the actual extinction could be underestimated. 

From Balmer decrement, the line ratio Ha/Hfi, the derived 
color excess of stellar continuum is Es{B-V) = 0.31 15, which 
is twice as the one obtained from spectral fitting. This indi- 
cates differential attenuation between the dusty starburst re- 
gion and older populations. 

6.3. IRAS F 13308+5946: A Possible Evolutionary Link 
between (type I) ULIRGs and QSOs 

The spectral synthesis gives a 500 Myr starburst history, 
which implies a triggering mechanism took place during such 
epoch. Galaxy interactions and mergers have timescales 
of a few xlO'^yr (Binney & Tremaine 1987). The similar 
timescales imply that the starburst was possibly triggered by 
galaxies merger. 

In §4 and §5, we demonstrated that IRAS Fl 3308+5946 has 
possibly experienced an IR QSO phase since ^^300 Myr ago. 
Foflow Hao05, we plot IRAS Fl 3308+5946 on the diagram 
which shows the relation between IR luminosity and the bolo- 
metric luminosity measurement (represented by ALa(5100); 
Kaspi et al. 2000) in Fig. E IRAS Fl 3308+5946 locates 
at a transitional position between IR QSOs and PG QSOs. 
The tight correlation followed by PG QSOs and NLSls sug- 
gests their IR luminosities are associated with the optical 



through central AGNs, while IR QSOs have IR excess pow- 
ered by other mechanisms rather than AGN. As for IRAS 
Fl 3308+5946, we attribute the IR excess to the ongoing star 
formation. In Fig. [TOl which shows the IR spectral index 



q;(60,25) (defined as a(Ai,A2) 



log[F(A2)/F(Ai)] 
log(A2/Ai) 



) versus the 



IR excess Lir/Lsooo, it also locates between IR QSOs and PG 
QSOs. The IR excess increase as the dust temperature de- 
crease, suggesting that star formation activities become more 
and more important in powering IR luminosity of IR QSOs 
and ULIRGs, because the temperature of dust heated by stars 
is lower than that by AGN. 

On the other hand, Hao05 gives statistical median Mbh val- 
ues of IR QSOs and PG QSOs, which are 4.9xlO^M0 and 
2.1x10*^^^0. Thus, IRAS Fl 3308+5946 has both BH mass 
and Eddington ratio similar to those of PG QSOs. Further- 
more, as a type 1.5 galaxy, the featureless continuum and 
emission-lines from the AGN may be partially attenuated by 
the dusty torus, and the AGN itself may be powering a classi- 
cal QSO, for the galaxy is an /-band quasar. 

Therefore, IRAS Fl 3308+5946 may evolve into an optical 
QSO when the starburst ceases and the nuclear dust is dis- 
sipated by the radiative pressure of AGN. In Fig.|9]and[T0l 
the positional arrangement of ULIRGs, IR QSOs and classi- 
cal QSOs implies an evolutionary sequence: ULIRGs — > type 
I ULIRGs / IR QSOs -^ (IRAS Fl 3308+5946) -^ PG QSOs. 
As the //-band luminosity is ^4L*, it support the evolutionary 
scenario proposed by Colina et al. (2001) that ULIRGs gen- 
erated by two or more massive (> L*) galaxies would evolve 
into QSOs (also supported by Sanders et al. 1988 and Lutz et 
al. 1999). 



IRASF13308+ag6» 



log XL,(5100A)/Ls 

Fig. 9. — IR luminosity vs. the bolometric luminosity measurement, 
AL;\5100A. The red open circles represent IR QSOs; the green filled circles 
represent PG QSOs; the brown filled triangles represent NLSls; the solid 
line represents the linear regression for all PG QSOs and NLSls. These data 
points are reproduced with the data from HaoOS. IRAS F13308+5946 is rep- 
resented by the blue star. 



Finally, through r-band surface photometry, we find the out- 
skirt surface brightness profile follows the de Vaucouleurs law 
and the central region follows a point spread function (PSF). It 
is consistent with other studies which show that some (maybe 
even higher fraction) ULIRGs may undergo QSOs phase in 
their evolutionary history before they settle down as ellipticals 
(Zheng et al. 1999; Arribas et al. 2004; Dasyra et al 2006b). 

7. SUMMARY 
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IRAS FY330ei-^946 



1 10 100 1000 10000 

^ir/^sioo 

Fig. 10. — IR spectral index o(60,25) vs. the IR excess, Lir/Ljioq. The 
purple open triangles represent 15 ULlRGs, which are selected from Kim et 
al. (1998b) and have SDSS spectra; other data points are also from Hao05 and 
have the same meaning as Fig.|9l the three rec(open circles enclosed by green 
diamonds represent IR QSOs which are also PG QSOs. IRAS F13308+5946 
is represented by the blue star and marked with a lower Umit. 



We carry out a study based on stellar population synthesis 
result of a type 1 LIRG, IRAS Fl 3308+5946. We find that: 

1. The cross-identification immlRAS, 2MASS, SDSS and 
FIRST, combining with the radio-FlR correlation, confirms 
that the LIRG IRAS Fl 3308+5946 is the optical quasar from 
SDSS observation. 

2. With sub-solar metallicity Z=0.008 and power-law in- 
dex a = -2.0, stellar population synthesis shows that the host 
galaxy has a recent starburst history since 500 Myr ago. 

3. We estimate the past IR luminosity during the starburst 
epoch. We find it has probably experienced a type 1 ULIRG 
phase from ^ 300 Myr ago. When the star formation activity 
weakened recently, the IR luminosity decrease to the present 



level (1O"'^*'L0) as a LIRG. Nuclear starburst and AGN ac- 
tivity both contribute to the IR luminosity budget, with ^-^70% 
from starburst. 

4. The SMBH mass Mbh = 1.8 x IO'^Mq and the Edding- 
ton ratio Lboi/i'Edd = 0.12 are both consistent with PG QSOs, 
suggesting a potential classical QSO. 

5. It locates at the transitional position between IR QSOs 
and PG QSOs on the plots Lir versus ALa5100 and a(60,25) 
versus Lir/Lsioo. Combining the results above, we conclude 
that IRAS Fl 3308+5946 is probable an evolutionary ti-ansi- 
tion object from a type 1 ULIRG to a PG QSO. 
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